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ABSTRACT 


We compare the performance of mass estimators for elliptical galaxies that rely on 
the directly observable surface brightness and velocity dispersion profiles, without in¬ 
voking computationally expensive detailed modeling. These methods recover the mass 
at a specific radius where the mass esti mate is expected to be least sensitive to the 
anisotropy of stellar orbits. One method (IWolf et al.ll2010ll uses the total luminosity- 
weighted velocity dispersion and evaluates the mass at a 3D half-li ght radius r 1 / 2 , i.e., 
it depends on the global galaxy properties. Another approach (IChurazov et al.ll2010tl 
estimates the mass from the velocity dispersion at a radius i ?2 where the surface 
brightness declines as i?~ 2 , i.e., it depends on the local properties. We evaluate the ac¬ 
curacy of the two methods for analytical models, simulated galaxies and real elliptical 
galaxies that have already been modeled by the Schwarzschild’s orbit-superposition 
technique. Both estimators recover an almost unbiased circular speed estimate with 
a modest RMS scatter (< 10%). Tests on analytical models and simulated galaxies 
indicate that the local estimator has a smaller RMS scatter than the global one. We 
show by examination of simulated galaxies that the projected velocity dispersion at 
i ? 2 could serve as a good proxy for the virial galaxy mass. For simulated galaxies the 
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total halo mass scales with a p {R/) as M v . 
RMS scatter « 40%. 
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1 INTRODUCTION 

Galaxy masses play a key role in our understanding of 
their formation and evolution. While observations of disc 
rotation curves allow the determination of spiral galaxy 
masses directly, the situation with early-type galaxies is 
more complex due to the lack of mass tracers on known 
orbits. Approaches of different levels of sophistication and 
generality have been developed for mass determination of 
elliptical galaxies. Dynamical modeling using the orbit- 
superposition method is considered to be the state-of-the- 
art technique and allows one to recover both the galaxy 
gravitational potential and the orbital structure with an 


accur acy of < 15% llThomas et al.l 120051 : iKrainovic et al.l 
l2005ll . Schwarzschild modeling is now widely used for mass 
measurements of supermassive black holes, for determina¬ 
tion of the total mass profile and its decomposition into 
luminous and dark matter compone nts as well as for con- 


straining the orbital structure fe.g., Gebhardt et al.1 

2003; 

CaDDellari et al. 

200t 

: Thomas et al. 2007bl. 20091. 

2011: 

McConnell et al 

2012 

2013: Rush et al. 201311. As such an 


approach only makes sense with high quality observational 
data allowing the determination of the high order line- 
of-sight velocity moments (namely, in addition to velocity 
and projected velocity dispersion also the third and the 
fourth order Gauss-Hermite moments fe.g.. [Gerhardl[l993l : 
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Ivan der Marel and Franxlll993lf or better yet the complete 
line-of-sight velocity distribution), it can only be applied to 
nearby galaxies. Moreover, numerical experiments show that 
due to intrinsic degeneracies not all the model parameters 
can be uniquely constrained even fr om the best available 
integral-field stellar kinematics (e.g. , iThomas et al.l l2007al : 
Ivan den Bosch and van de Venll2009f) . 

Large observational surveys at lower resolution have be¬ 
come a major tool for galaxy studies as they facilitate a num¬ 
ber of statistical investigations of galaxy properties. Here, 
exact galaxy mass determinations at different redshifts are 
critical for galaxy formation studies and for tracing the as¬ 
sembly of galaxy mass over time. With low-resolution infor¬ 
mation about galaxy photometry and kinematics, the usage 
of detailed dynamical modeling is impractical/not justified. 
It is desirable to have simple and robust techniques based 
on the most easily accessible observables that provide mass 
estimates with a modest and known scatter. 

Recently t wo simple mass estim ati on methods have 
been suggested dChurazov et al.ll20ld and lWolf et al J [ 2 01 (j| j 
which evaluate masses at a special radius from the surface 
brightness and projected velocity dispersion profiles without 
detailed modeling. Although these approaches recover the 
mass at some particular radius only and not the mass dis¬ 
tribution within that radius, these estimates could be used 
(i) to determine the non-thermal contribution to the total 
gas pressure when compared with the X-ray mass estimate 
at the same radius; (ii) to derive a slope of the mass profile 
when combined with the mass estimate from strong lensing; 
(iii) as a virial mass proxy. 

Here we compare the performance of simple mass es¬ 
timators on analytical models, simulated galaxies, and also 
on a set of real elliptical galaxies with high-quality kinemat- 
ical data and Schwarzschild modeling results. The paper is 
organized as follows. In Section 0 we provide a brief de¬ 
scription of the simple mass estimators. We present results 
of the tests on analytical models and a sample of simulated 
galaxies in Section 0 and on real elliptical galaxies in Sec¬ 
tion U The possibility of using these estimates as a proxy 
for the virial mass is discussed in Section [5] A summary of 
the accuracy and biases of the methods and conclusions are 
given in Section [SI 


2 MASS APPROXIMATION FORMULAE 


High resolution observations of the stellar surface brightness 
and line-of-sight velocity dispersion profiles of galaxies are 
the basis of various dynamical modeling techniques of dif¬ 
ferent levels of sophistication and generality. Unfortunately, 
all dynamical methods probably suffer from a common prob¬ 
lem: even with the best available kinematics and photometry 
they are usually unable to unique ly reso lve several degenera¬ 
cies (e.g., iDeionghe and Merritu 1992T I. The key problems 
are (i) the degeneracy between the total mass pr ofile and 
velocity anisotropy /3 aBinnev and Tremaine! l2008l l and (ii) 
the deprojection of the surface brightness profile I(R) into 
a three-dimensional lu minosity density for axisymmetric or 
triax ial systems (e.g., iRvbickil 1 19871 : I Gerhard and Binnevl 
1 199611 . 


The (scalar) virial theorem provides a straightforward 
way to estimate the mass (or circular velocity) of an ellipti¬ 
cal galaxy. For an isolated stationary spherical system in a 
logarithmic gravitational potential <f>(r) = V 2 ln(r) + const 
the circular velocity V c is related to the average (luminosity- 
weighted) line-of-sight velocity dispersion as 


Vc 


3(a 2 ) = 3 


f 0 °° al(R)I(R)RdR 
/” I(R)RdR 


( 1 ) 


The mass is then equal to M{< r ) = rVjf/G, where G 
is the gravitational constant. Under the above ideal assump¬ 
tions, this approach is rigorously independent of the velocity 
anisotropy. However, in practice, observed galaxies are not 
guaranteed (i) to have spectroscopic data over the full ex¬ 
tent of the systems, (ii) to be spherically symmetric and/or 
(iii) to have a stellar and dark matter distribution that sum 
up to an exactly logarithmic gravitational potential. 

To weaken some of these assumptions it is common to 
use the spherical Jeans equation which relates the veloc¬ 
ity anisotropy parameter /3(r), the luminosity density of the 
stars j(r) and their radial velocity dispersion ov(r): 


a I ■ 2\ . c,P ■ 2 
-(ja r )+2-ja r 



( 2 ) 


where /3(r) =1 — <rf/o/? (erg(r) is the tangential velocity 
dispersion; and in spherical symmetry ag = 0 $). 

In principle, one can derive M(< r) from the Jeans 
equation linking the 3D quantities j(r) and o>(r) to the 
observable surface brightness I(R) and projected velocity 
dispersion <J P (R) via 



ol{R)I(R) = 2 


ja 2 rdr 
\Jr 2 — R? 


( 3 ) 

( 4 ) 


for any given anisotropy profile. Unfortunately observational 
data alone do not allow one to constrain /3(r) without invok¬ 
ing sophisticated detailed modeling and additional observa¬ 
tional constraints. From considering relation © one might 
expect to find some characteristic radius R c har where the 
uncertainty in the circular speed estimate arising from the 
unknown anisotropy is minimal, thus mitigating the mass- 
anisotropy degeneracy at the expense of losing spatial reso- 
lution/in formation. The existence o f such a radius was first 
noted in iRichstone and Tremainel dl984lf . In this case the 
mass estimator can be expressed in the following form: 

Vc(Rchar) = fcffp, (5) 


where a p is some measure of the velocity dispersion designed 
so that the dependence of k on the orbital structure is as 
small as possible. 

This idea is the basis of the simp le d y namical ma s s scal - 
ing relations di s cussed in detail by IChurazov et al.l (l2010h 
and IWolf et ~ahl (l2010l i. In the remainder of the paper we 
concentrate on these two mass estimators and test their per¬ 
formance. 
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2.1 Local estimator 


The Churazov et al. estimator is derived from the station¬ 
ary non-streaming spherical Jeans equation under the as¬ 
sumption of a logarithmic gravitational potential <f>(r) = 
Vf ln(r) + const. In this case projected velocity dispersion 
profiles can be analytically derived for isotropic (/? = 0 ), 
radial (/3 = 1) and circular (/3 —> —oo) stellar orbits: 


crp°{R) = V c (r) 


VT+~o~+7 


< rc (i?) = V e {r) 


a™ d (R) = V c (r) 


2(1 + a + 7) 

1 


(6) 


\/(« + 7 ) 2 


+ <5-1 


where 

din/ _ diner 2 d 2 \n[Ia p ] 

a= ~d\nR’ 7 = _ dlnR ’ = d(\nR) 2 ' ^ 

These equations are only exact for a logarithmic po¬ 
tential, that is, if the circular speed is independent of ra¬ 
dius. However, on scales where the circular speed varies suf¬ 
ficiently slowly the projected velocity dispersions for /3 = 
0 , — 00 ,1 approximately follow relations ©> which allows 
one to invert eq. © and form estimates of V c for given 
anisotropy values (V), 130 , T4 clrc and V r c rad ) via local proper¬ 
ties of I(R) and <j p (R). 

For typical galaxy potentials, t here exists a character - 
istic radius R c har (called R a we et by IChurazov et al. H), 
where the circular speed estimate does not strongly depend 
on /3. Ideally, this radius would be defined as the point of in¬ 
tersection of V). 130 , V(, circ and F ( rad , i.e., where circular speed 
estimates for different /3 give exactly the same values. In real 
galaxies I 4 130 , VC' and V™ d are not guaranteed to inter¬ 
sect at one radius and _R awee t is defined as the radius where 
these three curves are most close to each other, i.e., where 
the circular speed estimate is largely independent of /3. In 
this paper we calculate 7/ awe et as the radius at which the 
standard deviation among the three curves is minimized. As 
seen from the equations above, Rawest depends only on the 
local properties of the observed surface brightness I(R) and 
projected velocity dispersion o- p (R) profiles. 

For massive elliptical galaxies the spatial variation of 
a p (R ) is typically much smaller than that of /(/?), i.e., 7 <C 
a and <5 -C a. Hence, the estimates of circular speed derived 
from (inverted) eq. © can be approximated by: 

V r “ 0l “ = a p (R)C^TI 


where the surface brightness declines as _R~ 2 . For general 
spherical models R s weet is expected to lie not far from R 2 
which, in turn, is close to the projected half-light radius 
(sometimes called the effective radius) Ri/%. Thus, in prac¬ 
tice, one can use R 2 instead of R s we et as the characteristic 
radius for the circular speed estimate. In this case the sweet 
radius is derived from the logarithmic slope of the surface 
brightness profile only. So we may consider the following 
three local circular speed estimators 

• estimator LI: Vjf^-Rsweet) where I 4 130 is given by the 
first of eq. © and .Rsweet is the radius at which the three 
estimators of V c in eq. © are closest to one another. 

• estimator L2: V)( so (/? 2 ) where Vj , 130 is the same as in 
LI and R 2 is the radius at which the logarithmic slope of 
the surface brightness distribution is —2 (a = 2 in eq. ©)■ 

• estimator L3: Vj. 130 ’ 3 (/? 2 ) where V ). 130 ’ 3 j s given by the 
first of eq. © and R 2 is the same as in L2. 


2.2 Global estimator 

The mass estimator suggested by IWolf et ~akl (|2010| J is de¬ 
rived from the spherical Jeans equation combined with the 
scalar virial theorem under the assumption of a fairly flat 
line-of-sight velocity dispersion profile a p (R). Wolf et al. 
show that the mass uncertainty arising from the unknown 
anisotropy of stellar orbits /3 is minimized at the radius r 3 
where the 3D luminosity density profile j(r) decays as r~ 3 . 
Within this radius the galaxy mass can be inferred from 
the luminosity-weighted projected velocity dispersion pro¬ 
file (a 2 ), averaged from 0 to 00 : 

M(< r 3 ) ~ 3 G~ 1 (< 7 2 )r 3 , (9) 

or in terms of the circular speed (hereafter, the global esti¬ 
mator Gl) 

Vf (r 3 ) ~ 3(c 2 ). (10) 

As seen from these equations, the mass estimator depends 
on the global properties of the observed surface brightness 
I(R) and the projec ted velocity disper sion cr p (J?) profiles. 

As discussed in IWolf et al] (I 2010 I) . for a wide range of 
stellar distributions (e.g. for exponential, Gaussian, King, 
Sersic profiles) r 3 is close to the 3D half-light radius 77 / 2 , 
which is in turn approximately equal to |i?i/ 2 , where R 1/2 
is the projected half-light radius. So the circular speed esti¬ 
mate can also be expressed as 

Vo (ri/ a ) ^ 3 (<t 2 ) (11) 

(estimator G2) or 


VC’ s = a p (R)J 2^ + i (8) 

V ot 

K rad ’ 3 = a P (R)V^C I. 

Note that for a = 2 the relation between the circular 
speed and the observed velocity dispersion is the same for 
isotropic, circular, and radial orbits and in fact this equality 
exten ds to all systems with constant anisotropy j3 llGerhardl 
Il993l h i.e., for the ideal case of a nearly flat projected veloc¬ 
ity dispersion profile, the characteristic radius Rsweet = R 2 


Vc (f fli/a) - 3 Cl) (12) 

(estimator G3). 

These mass estimators are intended for use with galax¬ 
ies that (i) are spherical, (ii) are non-rotating, (iii) have spa¬ 
tially resolved kinematics over the entire galaxy and (iv) 
have a projected velocity dispersion that does not vary with 
radius. 

Both the Churazov et al. and the Wolf et al. mass esti¬ 
mators are based on the spherical Jeans equation and only 
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Table 1. Main properties of the Churazov et al. (local) and the Wolf et al. (global) estimators. 


Estimator 

Assumptions 

Data 

Formulae V c ?{R c har) = 

Local 

(Ll, L2, L3) 

dynamical equilibrium 
spherical symmetry 

no rotation 

logarithmic potential 

log-slope of I(R), 
G'p(Rsweet') Of CTp(R 2 ), 
log-slope of <j p (R ), 
d 2 ln[/cTp]/d(ln R ) 2 

Ll: V)) ( fisweet ) — [1 + 0 :( /isweet. ) T 7( Rsweet )] ' Op (Rsweet ), 

where R SW eef. F c iso « F c circ « F c rad (eq. 10). 

L2: 2 (R 2 ) = [3 + 7(R2)] • 

where it 2 : I(R) oc R~ 2 . 

L3: V 2 (R 2 ) = 3 -o 2 (R 2 ). 

Global 

(Gl, G2, G3) 

dynamical equilibrium 
spherical symmetry 

no rotation 

flat cr p (R) 

<jp(R) over entire galaxy, 
deprojection of I(R) 
or determination of R1/2 

Gl: F c 2 (r 3 ) = 3© 2 ), 

where 7*3: j(r) oc r ~ 3 . 

G2: y c 2 (r 1/2 ) = 3( C r 2 ), 

where rq /2 i s the 3D half-light radius. 

G3: V c 2 (| 7 * 1/2 ) = 3(cr 2 ,) 

where R\/2 is the projected half-light radius. 


allow the determination of the total mass within a sphere of 
some particular characteristic radius Rchar, at which the 
mass estimate is not sensitive to the anisotropy. Despite 
these similarities, the final equations and expressions for 
Rchar look quite different. The crucial difference between 
these mass estimators is the following: the Churazov et al. 
14 -estimate depends only on local properties of the observed 
profiles (a p (R) and the log-slopes of I(R) and a p (R)), while 
the Wolf et al. formula requires averaging of the velocity 
dispersion over the entire extent of the galaxy and the de¬ 
termination of the projected half-light radius, i.e., it depends 
on the global galaxy properties. The main properties of the 
Churazov et al. and the Wolf et al. estimators are summa¬ 
rized in Table [T] 


3 PERFORMANCE TESTS 

In order to compare the performance of the local and global 
approaches, we test them on analytical models, simulated 
galaxies and real elliptical galaxies, the latter having also 
been analyzed using Schwarzschild modeling. 


3.1 Analytical models 

To set up the analytical models we solve the Jeans equa¬ 
tion numerically for a set of assumed 3D luminosity den¬ 
sity, anisotropy and circular speed profiles. The 3D lumi¬ 
nosity density comes from the spherical deprojection of the 
Sersic profile I(R) = /(i?i/ 2 )exp [—&„ (^R/R 1 / 2 ) 1/n — ljj 
dSersidn~968l J. where b n obeys T(2 n) = 27(2 n,b n ). T(2n) is 
the complete gamma function and 7(2 n, b n ) is the lower in¬ 
complete gamma function. We assume an anisotropy /3(r) = 
r c + /3iTa) / ( r c + r £). Here c is a concentration param¬ 
eter and r a is some characteristic anisotropy radius, where 
the anisotropy changes from its central value, /3 1 , to its outer 
value, /3 2 . In real early-type galaxies the distribution of stel¬ 
lar orbits is believed to be mildly tangential at the very 


center dThomas et al.ll2014l ). close to isotropic at small and 
moderate distanc es from the center and more radial in th e 
outer parts (e.g.. [Gerhard et al.ll200il : iThomas et al. IB)- 
Thus we vary /3i from —0.7 to 0.0 and /3 2 from 0.1 to 0.5. 
We vary r a and c in the range 0.1J?i / 2 ^ r a < 3.2 R 1/2 and 
0.1 < c ^ 5.1, respectively. 

The circular speed is parametrized as 


Vc {r ) oc Pi 2 


r 2 + r£ 


+ V2 


2 ln( l +r/r s ) -r/(r + r s ) 
r/r s 


(13) 


to mimic mea s ured 14-profiles (e.g., I Gerhard et al.l 1200 ll ; 
iBreddels et al.l 120131 ). In other words, the gravitational 
potential of a model spherical galaxy is assumed to 
be a combination o f the cored logarithmic potential 
dBinnev and Tremaine! 120081 ) and the potential which is 
si milar in shape to a Nav arro-Frenk-White potential 
dNavarro, Frenk fe Whitelll995l ). Let us note that this is a 
simple analytical representation of the circular velocity and 
the free parameters are chosen to cover the wide range of 
observed shapes and to probe different parts of real circular 
speed profiles (rising, roughly flat and decreasing 14). The 
parameter r s varies from 0 to 90 half-light radii, r c changes 
from 0 to 12 Ri/ 2 - 

The luminosity-weighted average of the projected 
velocity dispersion profile (cr 2 ) is calculated over 
[0.1 Ri /2 ; 10Ri / 2 ]. Figure □ illustrates all steps of the 
analysis and shows typical profiles considered for analytical 
models. 


3.1.1 Ideal models 

First we apply the estimators described in Section [2] and 
Table [T] to two ‘ideal’ galaxies, which meet all the assump¬ 
tions used to derive the local and global formulae: dynami¬ 
cal equilibrium, spherical symmetry, no streaming motions, 
constant anisotropy, flat V c (r) for the Churazov et al. for¬ 
mula or constant a P (R) for the Wolf et al. estimator. We call 
these two models the ‘ideal local’ and ‘ideal global’ galax¬ 
ies, respectively. Figure [2] shows typical profiles, derived for 
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0.1 1 10 


r /^l/2> R/^l/2 



0.1 1 R 2 r 3 10 


r /^l/2’ ^/^l/2 


Figure 1. Profiles for a typical analytical model (here, n = 4 Sersic model). Left: assumed surface brightness I(R ), anisotropy /3(r) and 
circular speed V c (r ) as functions of R/R 1/2 (or r/ R.] / 2 )- where R \/2 is the projected half-light radius. Right: profiles used to get the 
simple 14-estimates. Logarithmic slopes of j (r) (obtained from the spherical deprojection of the surface brightness) and I(R) are shown 
in the upper panel as dark green and purple curves respectively, r 3 , R Bwee t and R '2 are marked as dark green dashed, blue dashed and 
purple dotted lines. The velocity dispersion profile cr p (R) from the spherical Jeans equation is shown as a grey line in the lower panel. 
The simple circular speed estimates are shown as open and filled symbols of different colors. The filled dark green square is for the global 
Vf-estimate at r 3 (Gl). In this model r ^/2 is within ~ 1% of ^Ri /2 so we plot G2 and G3 together as a filled cyan square. The open 
blue square shows the local estimate at f? SW eet (LI) and the open purple square at R 2 (L 2 ), the open purple star shows the simplified 
version of the local estimator (L3). 


an n = 4 Sersic model with constant anisotropy /3 = —0.5, 
along with the derived local and global 14-estimates. The 
open blue square shows the local estimate at JJsweet (LI), 
the open purple square at R 2 (L2) and the open purple star 
is the simplified version of the local estimator (L3). Simple 
global estimates are shown by filled symbols: the dark green 
square is for the global 14-estimate at r 3 (Gl), the filled cyan 
square at rj/ 2 ~ f ^ 1/2 (G2/G3). As expected both estima¬ 
tors work well when applied to ‘ideal’ galaxies which meet 
all the corresponding assumptions. In the ideal local case, 
the global Wolf et. al estimators slightly underestimate the 
true circular speed with typical deviation of ~ —3%. Like¬ 
wise, in the ideal global case, the local estimators tend to 
overestimate the true 14 by ~ 3% corresponding to mass 
errors of =f6 %. 

3.1.2 Grid of analytical models 

We explore ~ 30000 analytical models, described by the 
Sersic surface brightness profile with index 2 < n < 20, 
anisotropy profile increasing slowly with radius and circu¬ 
lar speed characteristic for (i) dark-matter dominated dwarf 
spheroidal galaxies (14 growing with radius; ~ 50% of the 
models) and (ii) massive elliptical galaxies (14 roughly flat 
or decreasing slowly with radius; ~ 50% of the models). 
However, we do not aim to explore the whole parameter 


space: the idea is to understand how sensitive the estima¬ 
tors are to the assumption of a flat V c (r) or <r P {R) and to 
varying anisotropy. The resulting histograms for the local 
(upper row) and global (lower row) estimators are shown 
in Figure [3] The RMS scatter for the global estimators is 
almost twice as large as for the local one. 

As the Churazov et al. derivation assumes an logarith¬ 
mic gravitational potential (i.e., flat 14(f)) it is important to 
test how the accuracy of the estimator depends on the slope 
of the circular speed profile. Indeed, we find for the local 
estimator LI that there is a clear correlation between the 
fractional deviation A of the estimated circular speed from 
the true one and the logarithmic slope of the true V c (r) at 
Rsweet- A ~ k x (—din 14 2 /dlnr) (Figure [4] panel A). To 
understand this trend, let us consider for simplicity pure cir¬ 
cular orbits. If the circular speed increases (decreases) with 
radius, then stellar velocities at r > R are larger than at 
r = R, resulting in an increase (decrease) of a p compared to 
our reference case of a spatially constant 14- This, in turn, 
results in a systematic overestimation (underestimation) of 
the circular speed. The details of the mean correlation be¬ 
tween A and — din 14 2 /dln r depend on the sampling of the 
parameter space. However the main trend for the local esti¬ 
mator seems to be rather universal. For growing/decreasing 
14 near R S weet (which is close to R2) the method tends to 
overestimate/underestimate the true value of 14 by a factor 
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Sersic index n = 4, anisotropy (3= -0.5 Sersic index n = 4, anisotropy /3= —0.5 




r /R 1/2 . r /R, / 2 r / R i/ 2 > R / R l/2 

Figure 2. Circular speed estimates for an ‘ideal’ model galaxy with flat V c (r) (left side) or flat cr p (R) (right side). The surface brightness 
is described by the Sersic profile (Sersic index n = 4), and the anisotropy parameter is (3 = —0.5. The parameters n and (3 are chosen 
rather arbitrarily for demonstration purposes only. The upper panel shows the log-slopes of the 3D luminosity density 0 - 3 d (in dark 
green) and of the surface brightness a (in purple). The projected velocity dispersion, the true circular speed profile as well as simple 
Vc-estimates are shown in the lower panel. The symbols in the lower panels are the same as in Figure^ r%, R s weet and R 2 are marked 
as dark green dashed, blue dashed and purple dotted lines. 
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Figure 3. The histograms of the deviations of simple l/c-estimates from the true value V* rue for the local (upper row) and the global 
estimators (lower row) for model spherical galaxies (Section l3.1.2(l . The analytical models are described by Sersic surface brightness 
profile, anisotropy increasing slowly with radius, and circular speed profile that is similar in shape to observed circular velocity curves. 
Deviations are calculated as A = ( V c — Vc rue )/Vc rue , where the estimated V c and true V* rue are taken at the same radius. The 
RMS-scatter is calculated relative to A. 
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Y(^sweet) 



<*(r 3 ) a(r 3 )+y(r 3 ) 



Y( r 3) 


Figure 4. Dependence of the error in the 14 estimate on properties of the true V r (r) and the observable I(R) and a p (Ji) profiles for 
model spherical galaxies (Section l3. 1.211 . Panels A: Deviation A of the estimated V c from the true one as a function of the log-slope of the 
true circular speed — dlnV^ / d\nr taken at the characteristic radius (f? SW eet for the local estimator LI and r '3 for the global estimator 
Gl). Panels B: A as a function of the log-slope of the projected velocity dispersion 7 = —d\na p /d\nR. Panels C: A as a function of 
the log-slope of the surface brightness profile a = — d\nl/d In R. Panels D: A as a function of a + 7. The horizontal dashed line shows 
the deviation averaged over the whole sample of analytical models. The crosses show the average deviation in a chosen bin. Only 1000 
randomly chosen realizations are shown for clarity. Lower row: Color-coded deviation A as a function of a and 7. Left panels show the 
local estimator LI and right panels show the global estimator Gl. 


of ~ 1 + 0.1 (—din Vf /d\nr) . For flat 14 near Rsweet the 
local estimator is largely unbiased when averaged over the 
parameter space covered by our grid of analytical models. 
When deviations are plotted against the logarithmic slope 
of the projected velocity dispersion 7 = — din a p /d in R at 
Rsweet (Figure 2] panel B), a similar pattern is observed. 
If ap grows with radius in the vicinity of -Rsweet, then the 
local 14-estimate overestimates the true circular speed. For 
flat or moderately falling observed velocity dispersion pro¬ 
files the local estimator seems to recover an almost unbiased 
14-estimate. The observed trend in A ( 7 ) is partly compen¬ 
sated by the opposite trend in A(a) (Figure [I] panel C) 


to make the coefficient y 1 + a(R SW eet) +7 {Rsweet) relat¬ 
ing a p to the circular speed (first equation in eq. ©) be 
close to \/3. The deviation A is almost independent of a + 7 
(panel D of Figure|4|) for the local estimator, and the spread 
in (a + 7 ) at Rsweet is quite small (1.9 < a + 7 < 2.5). For 
the local estimators L2 and L3 correlations between A and 
—dh\Vf / dh\r are quantitatively the same as for LI. 

The global estimator shows a more complex dependence 
on the slope of the circular speed radial profile (Figure [4] 
panel A), giving a noticeable negative bias for a flat cir¬ 
cular speed (for flat 14 (r) the observed velocity dispersion 
might vary significantly with radius). As expected, the Wolf 
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8 Lyskova et al. 


et al. formula works best for roughly flat line-of-sight ve¬ 
locity dispersion profiles. The global 14-estimate appears to 
be too large for a p rapidly increasing with radius (7 < 0.3). 
Large negative deviations are present for models with er p (R) 
showing a bump, i.e., when 7 (R) changes sign. Similar de¬ 
pendencies are observed for the global estimators G2 and 
G3. 

Figure [4] indicates that the local estimator should be 
applied with caution to systems with increasing velocity 
dispersion profiles and/or to systems that are described by 
growing circular speed profiles in the vicinity of R 2 . As the 
global estimator relies on global properties of the galaxies, 
it works well if the velocity dispersion does not change sig¬ 
nificantly with radius over the whole extent of the system. 
Roughly speaking, the Wolf et al. formula is appropriat e for 
dwarf spheroidal galaxies (see iKowalczvk et all 120131 who 
have tested the global estimator on a sample of simulated 
dSph) and for a subset of elliptical galaxies with approxi¬ 
mately flat velocity dispersion profiles, while the local esti¬ 
mator works for elliptical galaxies in general. A similar con¬ 
clusion is reached when analyzing subsamples of analytical 
models with V c (r ) resembling that of (i) dwarf spheroidal 
galaxies and (ii) elliptical galaxies. 

It should also be noted that for large Sersic indices (n > 
8 — 10 ), typical for massive ellipticals sitting at the centers of 
groups or clusters, the log-slope of the surface brightness a is 
close to 2 over a wide range of radii, and in this radial range 
the true circular speed is well described by the isotropic one 
K 1S ° (eq. ©). 


3.2 Simulated galaxies 

From spherical models we now turn to a sample of ‘zoom-in’ 
cosm o logica l simulations of individual galaxies dOser et al] 
l 20 ld. 1201211 . which span a wide range in mass, 7 x 
10 11 Mq/i _1 < M v ir < 2.7 x 10 13 M q /i _1 , h = 0.72, where 
M v i r is the present day virial halo mass. Those simulated 
galaxies have properties very similar to observed nearby 


early-type galaxi es as studied by the ATLAS 3d project 
dNaab et al.l2013l f. The typical effective radii of these galax¬ 


ies are ~ 2 — 15 kpc. We define the effective radius to be 
i?i/ 2 , that is, the radius of a circle which contains half of 
the total stellar mass of a galaxy witho u t intr oducing any 
cut-off (in contrast to lOser et 2 ] | 2 Q 10 i . I 2012 II who define 
R \/2 as a radius which encloses half of the projected stellar 
mass within 10% of the virial radius). The typical axis ra¬ 
tios q (ellipticity e = 1 — q), calculated as the square root 
of eigenvalues of the diagonalised inertia tensor within the 
effective radius, are ~ 0.5 — 1. The anisotropy of the simu¬ 
lated galaxies is mildly tangential or close to isotropic at the 
center and be comes rad ially anisotropic with /3 ~ 0.2 — 0.4 
at large radii (|Wu et alj|2014h . The surface brightness pro¬ 
files of simulated galaxies are well discribed by cored Sers ic 
models with typical Sersic indices n > 10 (WoetaL 2014). 

_The local estimators have been tested in Lyskova et al.l 

(| 2012 |) and we follow the same analysis procedure. Briefly, 
we first exclude satellites from the galaxy image and cal¬ 
culate the radial profiles J(R), a P (R) and the true circu¬ 
lar speed Vj rue (r) = GM(< r)/r in a set of logarithmic 


concentric annuli/shells around the halo center. In practice, 
such averaging over concentric annuli could be applied to 
2D kinematic maps derived from integral field unit observa¬ 
tions. At .Raweet we take U c lao as an estimate of the circular 
speed and calculate the deviation from the true value at this 
radius A = (V / 80 - K true ) /R c true . 

We apply the estimators to a subsample of massive 
(a p (Ri/ 2 ) > 150 km s -1 ) simulated galaxies and exclude 
merging (~ 3% of the total) and oblate galaxies seen almost 
along the rotation axis (~ 15% of the total). 

Figure [5] shows the fraction of galaxies versus the devi¬ 
ation A = ( V c — Vj*™ 6 ) /14 true for the local estimators LI, 
L2, and L 3. The left and right p anels have already been pre¬ 
sented in I Lyskova et al.1 (l2012l l (their Figure 8). Note that 
LI and L2 estimators show very similar results: the average 
deviation over the sample is close to zero and the RMS scat¬ 
ter « 5 — 6 %. So the radius R 2 , which is uniquely determined 
from the slope of the surface brightness profile, can be used 
instead of R swee t, which depends also on the log-slope of the 
projected velocity dispersion profile — dlno p /d\n R and the 
second derivative d 2 \n[Ia 2 ]/ d{h\R) 2 . 

To test the robustness of the Wolf et al. estimator on 
the sample of simulated galaxies we need to compute (i) 
the luminosity-weighted projected velocity dispersion {a 2 ) 
over some radial range [Rmin] Rmax], (ii) the characteristic 
radii r 3 , 77 / 2 and 7 R 1/2 where the anisotropy is expected 
to play a minimal role in the circular speed determination. 
The 3D half-light radius ri /2 and projected half-light radius 
R\/2 are computed as the radius of the sphere/circle which 
contains half of the 3D/projected stellar mass, respectively. 
Since in real observations information on the total light is 
often not available or difficult to estimate, we also determine 
an ‘observational’ effective radius R e ff coming from a Sersic 
fit I(R) oc exp 6 n (R/R e ff) 1 ^ n ^ to the observed surface 
brightness over the same radial range [Rmin', Rmax] used for 
the calculation of (a 2 ). We choose Rmin = 3 r so ft, where 
Tsoft ~ 400 h -1 pc is the softening length of the simulations 
and 3r so f t is the maximum radius where profiles could be 
affected by the softening. Finally, R m ax = 0.1R„; r , where 
R v ir is the halo virial radius (= R 200 , the radius where the 
spherical overdensity drops below 200 times the critical den¬ 
sity of the Universe). 

Figure [6] shows the perfomance of the Wolf et al. esti¬ 
mator at different radii: (Gl) r 3 (upper left panel), where 
« 3 _d = — d\nj(r)/d\nr = 3, (G2) 3D half-light radius 77/2 
(lower left panel), (G3) |x ‘observational’ effective radius 
(upper right panel) and (G3') |x projected half-light radius 
R 1/2 (lower right panel). While the global estimates V c (77) 
and 14 ( 77 / 2 ) are almost unbiased (when averaged over the 
sample of simulated galaxies), the average deviation of the 
global V c (|R e ff) is « 3 — 4% high. The RMS scatter for all 
cases is equal to ~ 7 — 10%, i.e., significantly larger than for 
the local estimators. As mentioned above, observed profiles 
at R < Rmin = 3 r so f t could be affected by the softening, 
and for the analysis of the Wolf et al. estimator we consider 
only those simulated galaxies for which the characteristic 
radius is larger than 2 Rmin = Or soft- This selection crite¬ 
rion effectively keeps only the most massive galaxies with 
roughly logarithmic gravitational potential (see Section 3.2 
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Local estimator 



deviation from the true V e , % 


Figure 5. The histograms of the deviations of the local \4-estimators from the true circular speed for simulated galaxies (Section |3. 2 1) . 
Vc(Rsweet) (LI) and Vc(R. 2 ) (L2) are calculated from eq. J 6 } at Rsweet (radius where V^ lso ~ V^ irc ~ V^5 ad ) and R 2 (where I oc R~ 2 ) 
respectively. V£(R 2 ) (L3) comes from eq. ©■ 
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Figure 6. The histograms of the deviations of the global V c -estimators from the true circular speed for simulated galaxies (Section 13.21) . 
Deviations are calculated at the following characteristic radii: (i) 7-3 (upper left panel), where j(R) oc r~ 3 , (ii) a factor times the 
‘observational’ effective radius resulting from the Sersic fit to the surface brightness over [3r so ^ t ; O.l-R^r] (upper right panel), (iii) 3D 
half-light radius r-y /2 (lower left panel) and (iv) a factor ^ times the projected half-light radius R 1/2 (lower right panel). 
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in iLvskova et al.|[2012f) in the sample. For these galaxies the 
virial theorem states that 14 at any radius is well approx¬ 
imated by 3 (a p ) and the exact value of the characteristic 
radius is not important. 

If we vary Rmax and R' max within reasonable limits, the 
values of _R e ff and (cr p ) for individual galaxies do change, 
but not dramatically, and the average deviation A remains 
practically the same (|A| < 3 — 5%) with RMS scatter ~ 
8 - 10 %. 

Tests on model analytical galaxies fSection !3.1.211 show 
that for the local estimator the deviation A correlates with 
the logarithmic slope of the true circular speed. For the ris¬ 
ing I 4 (r) the local approach tends to overestimate the true 
circular velocity. For the sample of simulated galaxies we 
also observe such a trend in the local estimator LI (Fig¬ 
ure [7] Panels A), but not as strong as for the model spheri¬ 
cal galaxies. For the local estimators L2 and L3 the trend is 
similar. In the probed region of 7 the deviation A is fairly 
flat (Figure [7] Panel D) as is for the spherical analytical 
models (Figure H] Panel D). 

No clear correlation between simple local and global es¬ 
timates and the logarithmic slopes of the projected velocity 
dispersion 7 or the log-slope of the surface brightness a or 
the axis ratio is found. 

In Figure [ 8 ] we show the deviation of simple estimates 
from the true circular speed as a function of the virial 
halo mass M v i r . Dots represent the deviations for indi¬ 
vidual galaxies and lines are the least-square linear fits of 

A[%] = a ■ lg ^ l0 1 °A/ i .’/t ~ 1 ) + b' for the local es¬ 

timators LI and L2 and the global estimator G1 the trends 
with the virial mass are quite weak. The Wolf et al. formula 
has been applied only to simulated galaxies with the charac¬ 
teristic radius R c har ^ 2 Rmin = 6 r so f t , so the probed range 
of masses is different for different characteristic radii. 


4 COMPARISON OF SIMPLE MASS 
ESTIMATORS WITH A 
STATE-OF-THE-ART ANALYSIS 

We now proceed to test the simple mass estimators on real 
early-type galaxies which have high-quality resolved stellar 
kinematical data and which have been studied in detail, 
using Schwarzschild modeling. The number of such galax¬ 
ies is constantly increasing, but the extent of the kinematic 
data for the majority of them is limited to « R e ff (an ‘ob¬ 
servational’ effective radius based on the Sersic fit) only. 
For a successful mass determination with the simple esti¬ 
mators described here, however, it is desirable to have spa¬ 
tially resolved kinematics at least out to ~ 1.2 — 1.5R e ff- As 
our target sample we therefore use early-type galaxies from 
a Coma cl uster survey, modeled using the Schwarzschild 
method bv lThomas et all (l2007blh and, in addition, the gi- 
ant e lli ptical g alaxy M87 from the Virgo cluster modeled by 
I Murphy et al.l (1201 il l . For our purposes we choose only those 
galaxies from the Thomas et al. sample where kinematic 
measurements are available out to > 1.5R e ft and the galaxy 
is slowly rotating in the sense that a p (R e g) > V ro t{Reff), 


where a p and V ro t are the projected velocity dispersion and 
rotation velocity, both measured along the major axis. 

Galaxies of our target sample (7 Coma galaxies + M87) 
are listed in Table [7] We use the surface brightness and the 
root-mean-square velocity ( Vrms{R ) = \/& P (R) + V xot (R)) 
to recover the circular speed via eq. © and Gi. For each 
galaxy we compare simple circular speed estimators with 
the best-fit circular velocity. The luminosity-weighted ve¬ 
locity dispersion ( a p ) (which is needed for the Wolf et al. 
estimator) is averaged over the radial range where the ob¬ 
servational data are available. As some galaxies show a mod¬ 
erate rotation, we test also a slightly modified version of the 
Wolf et al. formula where ( a p ) is substituted by the average 
luminosity-weighted RMS velocity (Vjjms)- 

The results are presented in Figure [9] Panel A shows 
the rotational velocity V ro t (blue triangles) and projected 
velocity dispersion a p (red triangles) measurements along 
the major axis with errorbars (E ro t and E CT correspond¬ 
ingly). The RMS velocity Vrms = v P + V xot with ob¬ 
servational errors Ejjms = y/E^ ot + EJ is plotted as black 
triangles. The inter polated curves (the in terpolation proce¬ 
dure is deesribed in lChurazov et aHl20ld j for the kinematic 
data are shown as solid lines. Shaded regions show the mea¬ 
surement uncertainties. Panel B shows the best-fitting circu¬ 
lar speed vf chw from the Schwarzschild orbit-superposition 
modeling (black curve) with lcr-uncertainties (grey shaded 
region); plotted in blue is the isotropic circular speed calcu¬ 
lated from eq. © where cr P (R) is replaced with Vrms(R)- 
The blue shaded region reflects the uncertainties associated 
with VftMS-measurements. Simple 14-estimates are shown 
as squares of different colors: the blue open square corre¬ 
sponds to the local 14-estimate at R swee t (LI), the purple 
open square at R 2 (L2), global estimates (calculated as (a p ) 
and (Vr M s)) a t §R e ft (G3) are shown as the cyan and dark 
green filled squares. The effective radius R e a is shown with 
an arrow, and the radius R 2 as a purple dotted line. As 
the ‘true’ (coming from the Schwarzschild modeling) circu¬ 
lar speed and observed line-of-sight velocity dispersion (or 
RMS velocity) profiles are close to being flat, both Churazov 
et al. and Wolf et al. estimators are expected to recover the 
circular speed reasonably well. Note that galaxies in the cho¬ 
sen subsample do not match perfectly all the requirements 
for using the simple mass estimators. Namely, most of them 
are flattened and some are slowly rotating. Nevertheless, we 
see that the simple 14-estimates agree well with Vf chw , es¬ 
pecially for slowly rotating galaxies (GMP 0144, GMP 3510, 
GMP 5279), where simple estimates (local 14-estimates at 
Rsweet and R2) almost coincide with the V c Schw (r ) value. Ta- 
ble[3]summarizes the comparison of simple 14-estimates with 
the circular velocity from the Schwarzschild modeling. The 
estimators used are listed in Column (1), the mean deviation 
A from V c Schw in Column (2). RMS scatter around A (in 

ly N _, (A i - A ) 2 

Column (3)) is calculated as W 1 ---. Column (4) 

shows the uncertainty E 0 f, s associated with measurement er¬ 
rors. For each galaxy E 0 (, s is calculated as \J\ + a + q-E rms 
(eq. ©). The column (5) lists the average errors on the 
best-fit circular speed from the Schwarzschild analysis taken 
at the corresponding characteristic radius. For real ellipti- 
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Figure 7. Observed correlations for simulated galaxies (Section l3.2t) . Panels A: Deviation A as a function of the log-slope of the true 
circular speed — din V^/dln r taken at a characteristic radius (Rsweet for the local estimator LI and r% for the global estimator Gl). 
Panels B: A as a function of the log-slope of the projected velocity dispersion 7 = — d In a^/d In R. Panels C: A as a function of the 
log-slope of the surface brightness profile a = — dln//dln R. Panels D: A as a function of a + 7 . Panels E: A as a function of axis ratio. 
The horizontal dashed line shows the deviation averaged over the whole sample of simulated galaxies. The crosses show the average 
deviation in a chosen bin. Left panels show the local estimator LI and right panels show the global estimator Gl. 
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Figure 8. Deviation of the estimated circular speed from the true one as a function of the virial halo mass (extracted from simulations). 
The local formula (left side) recovers the true circular speed for a wide range of virial masses almost equally well. The global approach 
works better for the most massive galaxies in the sample (M v { r > 3 • 10 12 Mq/i _ 1 ), than for the less massive galaxies. One can notice 
that the global estimator at different characteristic radii probes a different range of virial masses. When using the global estimator at 
Rchar = {r" 3 , v\/ 2 , |R e ff} we retain in the sample only the galaxies which have R c har > 2 Rmin — Qr so ft- As a result the three panels 
for Rchar = { r ’ 3 5 r i /2 5 §^eff} probe different subsamples at the low-mass end. In addition, the characteristic radii 7*3 and 77 / 2 > which 
are expected to be close to each other for commonly used analytical stellar light distributions (e.g., for Sersic models) sometimes differ 
by a factor ~ 2 — 3 in the simulated galaxies. 
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Table 2. Sample of real e lliptical g alaxies an alyzed using the Schwarzschild modeling (7 Coma galaxies from iThomas et alj ll2007bf) 
and M8 7 (NGC4486) from [Murphy et al .11 201 71 ). The effective radi i and ellipticities of the Coma galaxies are taken from lThomas"et**all 
(2007.3), the effective radius and the ellipticity of M87 comes from iKormendv et al.l (l2000l i . 


name 

other common names 

R e ff, arcsec 

ellipticity e at R e ff 

NGC 4957 

GMP 0144 

18.4 

0.256 

NGC 4952 

GMP 0282 

14.1 

0.315 

NGC 4908 

GMP 2417 

7.1 

0.322 

NGC 4869 

GMP 3510 

7.6 

0.112 

IC 3947 

GMP 3958 

3.3 

0.323 

NGC 4827 

GMP 5279 

13.6 

0.205 

NGC 4807 

GMP 5975 

6.7 

0.170 

NGC 4486 

M87 

194.41 

0.218 


Table 3. Simple V^-estimates and V c from dynamical modeling. The columns are: (1) - the simple 14-estimator; (2) - Deviation of the 
estimated 14 from ySchw resu ]( j n g from dynamical modeling, averaged over the sample of 8 galaxies (7 Coma galaxies ± M87); (3) 
- RMS scatter around the average deviation; (4) - average observational error at the characteristic radius normalised to ySchw. f.yj _ 
average error on dynamically derived ySch-w , :| ) ;| le characteristic radius. Since the bias for the local estimators is mostly driven by a 
single galaxy IC 3947 (the most compact galaxy in the sample), we exclude this galaxy and provide the results of the analysis of the 
remaining galaxies in the parentheses. 


estimator 

(i) 

A = (V c /V c Schw - 1),% 
(2) 

RMS scatter, % 
(3) 

(X ob s)/V c Schw , % 

(4) 

(^Sch w )/V c Schw ,% 

(5) 

LI [Vc(-Rsweet)] 

5.2 (3.2) 

6.6 (3.1) 

±5.8 (±5.7) 

+ 8 ' 1 /— 4.8 

(+ 8 -°/_4 9 ) 

L2 [V c (R 2 )} 

3.3 (1.2) 

7.5 (4.9) 

±6.2 (±6.3) 

+ 9 ' 6 /- 6.4 

(+ 9 ' 7 /- 6 . 7 ) 

G3 [yj3(a%) at f R efI ] 

2.0 (1.1) 

7.4 (7.6) 

±5.7 (±5.8) 

+ 8 - 4 /- 5.4 

( + 8 ' 4 /- 5 . 6 ) 

\/^(^RMs') at t^eff 

8.9 (7.3) 

7.8 (6.9) 

±5.7 (±5.8) 

+ 8 - 4 /- 5.4 

( + 8 ' 4 /- 5 . 6 ) 


cal galaxies the characteristic radii ri /2 and r 3 may not be 
available, and their determination often requires additional 
assumptions, so in Table Owe provide results for the global 
estimate at |-R e ff only. Furthermore, elliptical galaxies in the 
target sample have roughly flat circular speed and projected 
velocity dispersion profiles, so the performance of the global 
estimator at all characteristic radii (r 3 , tt /2 and f-R e ff) is 
expected to be similar. The average deviation of the global 
Vdr?,) and V c {r 1 / 2 ) are ~ 2.6% and ~ 1.3% high, and the 
RMS deviations are ~ 6.0% and ~ 8.5%, correspondingly. 
Here, the radii r 3 and ri /2 are obtained from the spherical 
deprojection of the observed surface brightness profile and 
its Sersic fit, correspondingly. 

Figure [TT] shows the observed correlations for the target 
sample of real early-type galaxies. Again, there is no clear 
correlation between simple 14-estimates and axis ratios of 
galaxies, although the selection criterion on the rotational 
velocity ( a p > V ro t at the effective radius) leaves only rela¬ 
tively round galaxies. 

For galaxies with a significant amount of rotation (i.e., 
Vrot(Reff) > cr p (J?eff)) simple methods provide notably bi¬ 
ased 14-estimates. Figure [l^] illustrates the comparison of 
simple estimates with the Schwarzschild dynamical model¬ 
ing of four additional early-type galaxies from the Thomas 
et al. sample, which have kinematics along the major axis 
available out to > 1.5f? e ff and cr p (R e ff) < V ro t(Res)■ When 
the RMS velocity {Vrms = y/cr% + V POt ) is used for deriv¬ 


ing the circular speed of a galaxy, then the final 14-estimate 
overestimates the ‘true’ circular speed (coming from dy¬ 
namical modeling) by 15 — 30%. When one does not ac¬ 
count for rotation and uses only projected velocity disper¬ 
sion for estimating 14, then the final result is biased low by 
15 — 30%. The truth seems to lie somewhere in between. The 
14-estimate inferred from Vrms is expected to be almost 
unbiased onl y after averaging o ver galaxies with random 
inclinations ( Lyskova et al.ll2014l ). The investigated galax¬ 
ies with a significant amo unt of rotation are ma inly oblate 
galaxies seen edge-on (see IThomas et aD [20073) . And it is 
not surprising that for such objects with measured kinemat¬ 
ics along the major axis only the estimator in a form of 
V?(R char ) = with £ ~ 0.5 would give more 

sensible 14-estimate. As a crude approximation, we can as¬ 
sume that along the minor axis the rotation would vanish 
and the velocity dispersion measurements remain almost the 
same. Then averaging Vrms over major and minor axes 
would result in Vrms ~ \/<^p + 0.514 2 of where a p and V r ,,t 
are measured along the major axis. 


Rotation seems to be the main factor which drives the 
bias of simple circular speed esti mates. For this samp le of 
real early-type galaxies (11 from IThomas et al.l l200Tbl plus 
M87) the deviation (14 — Vf chw )/y^ ch,w correlates with the 
luminosity-weighted ratio 14ot/cr p measured within the ef¬ 
fective radius surprisingly well (Figure IT3ll . 
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Figure 9. Comparison of simple V>-es 1 i males with the circular speed (Vf chw ) coming from the dynamical modeling (Sect,ion |~Q). Panel 
A shows data on the projected velocity dispersion cr p and rotation velocity V ro t, both measured along the major axis, and RMS velocity 

(' Vrms = \J a p + V roii) as well as interpolated curves used to calculate the logarithmic derivatives. Panel B presents the circular speed 
resulting from the Schwarzschild modeling (black thick curve) with error bars (grey shaded region), the isotropic circular speed V" c lso (in 
blue) calculated via eq. ((6j) with ‘observational’ error bars determined from the measurement uncertainties on Vrms- The circular 
speed profiles for pure circular and pure radial orbits are shown as magenta and green curves respectively. Simple Vf-esti rnales are shown 
as squares of different colors: the blue open square corresponds to the local F c -estimate at f? swee t (LI), the purple open square at R 2 
(L2), global estimates (calculated as {a 2 ) and (Vf M g)) at ^R e s (G3) are shown as the cyan and dark green filled squares. The arrow 
and the dotted purple line show the effective radius i? e fj and R 2 respectively. 


5 MASS PROXY 

We also test on simulated galaxies whether the simple cir¬ 
cular speed estimates could be used as a proxy for the virial 
galaxy mass. The simple estimators rely on <t p (R, 2 ) and (cr 2 ). 
Let us plot these against the virial halo mass for each indi¬ 
vidual galaxy to see how well these quantities correlate with 
each other. 

Figure [14] shows the virial galaxy mass as a func¬ 
tion of the local value of the projected velocity dispersion 
R 2 (left side) and the luminosity-weighted average veloc¬ 
ity dispersion (of) 1 / 2 (right side). Colored squares depict 
the velocity dispersion and the virial mass for individual 


simulated galaxies (upper panels) and the straight line is 
the least-square log-linear fit to these data points. Devi¬ 
ations from the linear fit are shown in the lower panels. 
The virial mass M v i r (in Mgft -1 ) can be approximated 

by M V i r « 6.0 ■ 10 12 ( fJp( ' R ^ ) with RMS ss 38% or 
\200 kms -1 ) 

M vir « 4.5 • 10 12 [ . I with RMS « 37%. The 

dispersion measures cr p (R. 2 ) and (cr 2 ) predict halo mass with 
almost the same scatter, (a 2 ) is expected to serve as a 
good proxy for the virial halo mass for our sample of sim¬ 
ulated galaxies as the majority of them have almost loga- 
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Figure 10. (continue) 


rithmic gravitational potential which is well approximated 
by <h(r) « 3(cr 2 ) lnr + const at any radius according to the 
virial theorem. Curiously, the relation between supermassive 
black hole mass M bh and host-galaxy b ulge velocity disper¬ 
sion Mbh oc cr 4 ' 24 (lGiiltekin et al.ll2009l ) has approximately 
the same scatter (~ 31%) and the same exponent as the 
M v i r — o p (R, 2 ) and M v i r — (cr 2 ) relations. 

We also test whether CT p (i?i/2) and the luminosity- 
weighted projected velocity dispersion (cr 2 ) e measured 
within an aperture of radius R \/2 correlate with the virial 
mass. We find that the RMS scatter for cr p (Ri/ 2 ) is ~ 60% 
and for (<7 2 ) e the scatter is « 50%, i.e., noticeably larger 
than for <T p (i?2) and (<r 2 ). We have not experimented to de¬ 
termine whether <t p (J? 2) and (cr 2 ) are the best predictors of 
the halo mass. Therefore, we do not exclude that velocity 
dispersion measurement at/within other radii could work 
better. Note that in contrast to the estimators discussed 
in Sectional the halo mass estimators are empirical results 
only obtained for the sample of massive isolated simulated 
galaxies, with no dynamical justification. 


6 DISCUSSION AND CONCLUSIONS 

We have compared the performance of two simple and 
fast methods to evaluate masses of elliptical galaxies at a 
special radius where the mass estimate is largely insensi¬ 
tive to the anisotropy in velocity dispersion. Such methods 
could be useful for mass determination of large samples of 
galaxies with poor /noisy data when detailed investigation 
is not practical. A reliable mass estimate at a single ra¬ 
dius could also be used as an additional constraint for the 
Schwarzschild dynamical modeling thus reducing the range 
of gravitational potentials to be explored. 

One approach uses local properties of the galaxy - log¬ 
arithmic slopes (and sometimes curvature) of the surface 
brightness and velocity dispersion profiles and recovers the 
mass at a radius whe re the surface brightn ess declines as 
RT 2 (local estimator; IChurazov et al] 120101) . Another ap¬ 
proach uses the total luminosity-weighted velocity disper¬ 
sion (cr 2 ) and evaluates the mass at the radius where the 
3D luminosity density j(r) oc r -3 which can be related to 
3D and projected half-light radii as rz « rj/ 2 ~ \R\/i 
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Figure 11. Observed correlations for real galaxies (7 Coma galaxies + M87). Panels A: Deviation A as a function of the log-slope of 
the circular speed from dynamical modeling taken at a characteristic radius (R swee t for the local estimator LI and ^ Reff f° r the global 
estimator G3). Panel B: A as a function of the log-slope of the projected velocity dispersion 7 = — din d 2 /d In R. Panels C: A as a 
function of the log-slope of the surface brightness profile a = —dlnl/dlnR. Panels D: A as a function of a+ 7 . Panels E: A as a function 
of axis ratio 1 — e. The horizontal dashed line shows the deviation averaged over the whole sample of elliptical galaxies. The histogram 
shows the average deviation in a chosen bin. Left panels show the local estimator LI and right panels show the global estimator G3. 


Table 4. Comparing the performance of different estimators for different samples. Best-perfoming estimators are shown in bold face. 




Local estimator 

Global estimator 

Applicable to 


slowly rotating elliptical galaxies 

slowly rotating stellar systems 
with nearly flat d p (R) 

Data 


log-slope of I(R ) near Rchar-, 

Gpy.Rchar) 1 

log-slope of d v (R) near R r hnr and 
d 2 \n[I dp)/ d(\n R) 2 (for LI) 

deprojection of I(R) 
or determination of R 1 / 2 , 
dp(R) over entire galaxy 

(' vFiRchar ) = ko*) 

Rsweet Or R 2 

r 3 or r 1 / 2 or f-R 1/2 

dp 


dp (Rchar ) 


k 


1 T 0:( Rchar ) T 7 {Rchar) 
or 1 + a (# 2 ) 

3 

Estimator 


Ll [Vc (-Rsweet )] L2 [V c (R 2 )] 

G1 [V c (r 3 )] G2 p4(r 1/2 )] G3 [V c {\r 1/2 )\ 



spherical analytical models 

Average deviation 
RMS , % 

A, % 

1.78 ±0.02 2.50 ± 0.03 

3.66 4.48 

-4.00 ±0.03 -2.75 ± 0.04 -2.49 ± 0.04 

5.56 6.59 7.01 



simulated galaxies 

Average deviation 
RMS , % 

A, % 

0.0 ± 0.8 - 0.1 ± 0.8 

5.4 5.6 

0.7 ±1.3 0.1 ±1.6 3.4 ±1.7 

6.9 9.6 7.9 



real elliptical galaxies 

Average deviation 
RMS , % 

A, % 

5.2 ± 2.3 3.3 ± 2.7 

6.6 7.5 

2.6 ±2.1 1.3 ±3.0 2.0 ±2.6 

6.0 8.5 7.4 
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Figure 12. Compa r ison of simple Vfoestimates with the circular speed (V^ chw ) coming from the dynamical modeling for four galaxies 
from lThomas et al.1 (l2007lj ) which have kinematics available out to > 1.5it e ff and cr p (R e fj) < V ro t (R e ff)■ Panel A shows Vrms(R) = 
\JcTp ( Ji) + \\ ot ( fi) (in black), a v {R) (in red) and V ro t(R) (in blue) profiles measured along the major axis as well as interpolated curves 
used to calculate the logarithmic derivatives. Panel B presents the circular speed resulting from the Schwarzschild modeling (black thick 
curve) with error bars (grey shaded region), the isotropic circular speed calculated from the RMS velocity V c lso = VrmsV^ + a + 7 
(in blue, solid lines) and the isotropic V c which does not account for rotation V)!' so = a v V 1 a + 7 (in blue, dotted lines). The simple 
local Pc-estimates at H 2 (where a = din//din it = 2 ) are shown as open purple squares, the simple global estimates at ijit e fj - as filled 
squares. The arrow and the dotted purple line show the effective radius it e fj and R 2 respectively. The purple cross marks the average 
between simple local estimates. 


for a wide range of stellar distributions (global estimator; 
IWolf et al.1 1201(1 1. We test the accuracy and robustness of 
these simple mass estimators on analytical models, on a sam¬ 
ple of cosmological zoom-simulations of individual elliptical 
galaxies and on real elliptical galaxies that have already been 
analyzed by means of a Schwarzschild approach. 

We have found that: 

(i) for the analytical models both methods recover the 
true circular speed with high accuracy. For a grid of explored 
models the average deviation of the simple local W-estimate 
from the true circular speed is A ~ 2% with RMS scatter 
~ 4%. The global method gives A ~ —4% with ~ 6% scat¬ 
ter. Although the exact values of the average deviation and 


the RMS scatter depend on the sampling of the parameter 
space, the local estimator seems to be less sensitive to the 
assumptions under which it has been derived than the global 
one. 

(ii) We also examined massive ( <j v (R \/2 ) > 150 km s _1 ) 
simulated galaxies, excluding oblate objects seen almost 
along the rotation axis (~ 15% of the total). For these the 
local formulae LI, L2 recover an (almost) unbiased estimate 
of the circular speed with RMS scatter ~ 5 — 6%. The Wolf 
et al. relation also gives an almost unbiased measurement of 
V c at the radius where the 3D luminosity density declines as 
r~ 3 , with RMS « 7%. At |R e ft (where the effective radius 
R e ff is defined from the Sersic fit to the surface brightness 
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Local estimator LI [V c (R sweel )] Global estimator G3 [V c (4/3R eff )] 




Fig ure 13. Influence of rotation on the simple V^-estimates derived for the sample of real early-type galaxies (11 from [Thomas et al.l 
l2007bl plus M87). Deviation A as a function of V ro t/<?p , measured within R e g, is shown on the upper row. The lower row presents the 
color-coded A as a function of V ro t /and axis ratio. 




cr p (R a ), km s' 1 V<a^>. km s _1 

Figure 14. The virial mass of simulated galaxies as a function of the projected velocity dispersion at R 2 (left) and the luminosity-weighted 
RMS velocity dispersion (right). 


profile) the average global circular speed estimate is biased 
high by 3.4% and the RMS scatter is around 8%. For real 
elliptical galaxies _R e ff is subject to additional uncertainty 
(especially when the Sersic index is large) as its determina¬ 
tion depends on the radial range used for the analysis and 
applied methodology and this is another advantage of the 
local estimator. 

(iii) For a sample of eight real slowly-rotating ellipti¬ 
cal galaxies with a p (R e ft) < V r ot{Reff) analyzed with the 
Schwarzschild approach both methods show a remarkable 


agreement with the best-fit circular speed coming from the 
dynamical modeling. When averaged over the sample of 
eight galaxies our simple estimator overestimates the best- 
fit dynamical circular speed by 5%. This bias is mostly 
driven by a single galaxy (IC 3947) with the smallest R e g 
in the sample. When this galaxy is excluded from the sam¬ 
ple, the V r c (i? su ,eei)-estimator gives a circular speed estimate 
that is high on average by 3.2% relative to the vf chw with 
RMS scatter « 3.1%. The RMS scatter between our sim¬ 
ple estimates for the investigated sample of eight galaxies 
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is as 6.6% which is comparable to measurement uncertain¬ 
ties. The Wolf et al. estimator for the same sample gives a 
mean overestimate as 2% with slightly larger RMS scatter 
of w 7.4%. 

(iv) A galaxy appearing round on the sky could also 
be an intrinsically flattened system (e.g., an oblate galaxy 
viewed along the polar axis), in which case the simple esti¬ 
mates (and the Schwarzschild models) are expected to be¬ 
come less accurate. For samples of massive ellipticals, the 
contamination w ith flattened, face-on galaxies is expected 
to be small le.g.. [Emsellem at al.|[20llh . 

Tabic [3] provides the average deviation and the RMS 
scatter for different simple estimators resulting from the 
tests on spherical analytical models, simulated and real el¬ 
liptical galaxies. Note that the local estimators at _R swee t 
and show almost the same perfomance, suggesting that 
instead of searching for a radius -Rsweet by solving eq. m one 
can use the radius R 2 , where the surface brightness declines 
as R~ 2 , as the characteristic radius. 

The projected velocity dispersion value at the radius R 2 
where the surface brightness declines as R~ 2 seems to be 
a good proxy for the virial galaxy mass. M v i r (in Mq/i -1 ) 

can be approximated by M v i r as 6T0 12 ( — a p(^ 2 ^ [ w jth 

\200 kms -1 / 

RMS scatter as 40%. The scatter is comparable to the scatter 
observed when sj{a |) is used as a proxy for the virial halo 
mass. 

Wh ile we were writing this paper, we cam e across a 
paper of lAgnello, Evans and Romanowskvl ll2014h . who sug¬ 
gest, in particular, an ‘extension’ of the local approach for 
power-law total density profiles (ptot oc r~ a ) in the form 
of V 2 (Rm) = Ka 2 (Rcr), where R a is a radius where de¬ 
pendence on the exponent a is minimal, Rm is chosen in 
a way to minimize dependence of R& on the anisotropy 
profile and K is a dimensionless constant. For exam¬ 
ple, in case of the Sersic surface brightness with n = 
4 and the Osipkov-Merritt anisotropy profile the circu¬ 
lar speed of a galaxy can be estimated as V c (3AR e s) ~ 
1.67<7 p (1.1 5-R e ff). To apply the simple mass estima tor pro¬ 
posed by lAgnello, Evans and Romanowskvl J2014l 'l to our 
grid of analytical models and the sample of simulated galax¬ 
ies, we derived triplets (Ra, Rm, K) for the correspond¬ 
ing Sersic surface brightness fits and the Osipkov-Merritt 
anisotropy profiles f}(r) = r 2 /(r 2 + r 2 ) with r a = li? e ff, 3i? e ff 
and 10-Reff- For the sample of analytical models the average 
deviation is found to be ~ —7.4% and the RMS scatter 
is ~ 3.4%. The negative bias is also found for the investi¬ 
gated simulated galaxies: A = —10%, RMS ~ 5%. For real 
early-type galaxies the circular speed profile derived from 
the Schwarzschild modeling at large radii (R ~ 3J? e fi) be¬ 
comes quite uncertain making a comparison of the estimated 
Vc(Rm) at Rm > 3-Reff with the ‘true’ one not informa¬ 
tive. So from the tests on spherical analytical models and 
simulated galaxies we conclude that the simple estimator 
V 2 (Rm) = Ka 2 (R a ) seems to have no obvious advantages 
over the discussed global and local approaches. 
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